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ABSTRACT 

We present multi-frequency radio observations of the 2010 nova event in the symbiotic binary 
V407 Cygni, obtained with the Karl G. Jansky Very Large Array and spanning 1-45 GHz and 17-770 
days following discovery. This nova — the first ever detected in gamma rays — shows a radio light curve 
dominated by the wind of the Mira giant companion, rather than the nova ejecta themselves. The 
radio luminosity grew as the wind became increasingly ionized by the nova outburst, and faded as the 
wind was violently heated from within by the nova shock. This study marks the first time that this 
physical mechanism has been shown to dominate the radio light curve of an astrophysical transient. 
We do not observe a thermal signature from the nova ejecta or synchrotron emission from the shock, 
due to the fact that these components were hidden behind the absorbing screen of the Mira wind. 

We estimate a mass loss rate for the Mira wind of w 10~^ Mq yr"-'^. We also present the 
only radio detection of V407 Cyg before the 2010 nova, gleaned from unpublished 1993 archival VLA 
data, which shows that the radio luminosity of the Mira wind varies by a factor of >20 even in 
quiescence. Although V407 Cyg likely hosts a massive accreting white dwarf, making it a candidate 
progenitor system for a Type la supernova, the dense and radially continuous circumbinary material 
surrounding V407 Cyg is inconsistent with observational constraints on the environments of most 
Type la supernovae. 

Subject headings: Stars: individual (V407 Cygni) — binaries: symbiotic — novae, cataclysmic variables 
— Radio continuum: stars 
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Figure 1. Multi-frequency radio light curve for V407 Cyg. A variety of colors denote different frequencies spanning 1.4-45 GHz, as marked 
on the left side of the plot. The time of optical discovery of the nova event is 2010 March 10.8. 
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1. INTRODUCTION 

On 2010 March 10, the symbiotic star V407 Cygni 
hoste d the first nova explosion ever detec ted in gamma 
rays (jCheung et al.ll20li3t lAbdo et al.|[201Q) . This ex- 
plosion also constitutes the second example of an "em- 
bedded" nova thoroughly observed at a range of wave- 
lengths, presenting an opportunity to trace the inter- 
action of a blastwave with dense circumstellar material 
in detail (the first example was RS Ophiuch i, with re- 
cent outbursts in 1985 and 2006; Bode 1987; Evans et alj 
|2008[) . Ma ny parallels can be drawn between RS Op h and 
V407 Cyg (iMikolaiewskall2010l: iNelson et al.l[20T2h . and 
together they form a class of intriguing candidates for 
Type la super nova (SN) progenitors (Munari & RcnzinJ 
119921: lHachis"^l2003l : iPatat et aT]l2011al ). 

V407 Cyg is a D-type (dust-rich) symbiotic binary fea- 
turing a w hite dwarf (WD) an d a Mira giant star in a 
wide orbit ()Munari et al.lll99"ol). The giant pulsates with 
a period o f 763 days (jMeinungeii 119661 : iKolotilov et alJ 
119981 120031. the longest known period for a Mira symbi- 
otic system. There is some indication of a 43-year orbital 
perio d in the long-term optical light curve (jMunari et al.l 
|1990[ ). corresponding to a binary separation of 17 AU for 
a WD mass of 1 Mq and a giant mass of 1.5 Mq. 

V407 Cyg has undergone periods of increased activity 
throughout the past several decades, observed as peaks 
in the optical light curve and brightening e mission lines, 
particularly in th e years 1936 and 1998 (iHoffmeisteH 
1949; We ljnl [19731 IMunari et all [l990l: IKolotilov et al ' 
199ai2003n'Tatarnikova et al.ll2003bl:lEsipov et al.ll2012[ ) 



These outbursts lik ely mark an increase in the accretion 
rate onto the WD (IKolotilov et al.l [20031 ). and so do not 
qualify as bona fide novae, which are powered by a ther- 
monuclear runaway. 

However, the 2010 outburst was a true nova ex- 
plosion, displaying an unprecedented brightening of 
the optical light curve (reaching a B-band magnitude 
of ~8 mag, as opposed to 1 3.1 mag in the 1998 
outburst; Kolotilovelalj 120031 : iNishivama et al.l 120101 : 
IMunari et al.. .20111) and an expansion velocity of the 
ejecta of 3,200 km s"! (fB^lmnl: lAbdo et al.ll20Toh . 
The unique gamma-ray signature from this explosion 
has been attributed to particles accelerated by the in- 
teraction of the nova shock with the dense Mira wind 
(lAbdo et al.l[2OT0t IL ii et al.ll2011t IH ernanz fc Tatischefj 
I2012t iSitarek fc Bed narek 20127 !Martin fc Dubu£ll2012[ ). 
The optical prope rties of V407 Cyg's nova event have 
been describ ed by IMunari et al.l (|2011l ) and iShore et"ari 
(I20T1I [Ml. whil e the X - ray evolution wa s stud - 
ied bv IShore et al.l ([20111). lOrlando fc Drakj ([Ml) . 
iSchwarz et al.l ()201lD . and lNelson et al.l (|2012D . 

This paper presents complementary radio continuum 
observations of the 2010 outburst, obtained between 1- 
45 G Hz with the Karl G. Jansky Very Large Array 
(VLA: iPerlev et al.ll20Tll ). Before the 2010 outburst, no 
radio detections of V407 Cyg had ever been reported 
(although, see Section [7|), but shortly after the nova 
explosion V407 Cyg was strongly detected across the 



centimeter and millimeter bands ("Nestoras et al 



2010; 



2010; 



Gawronski et al.ll2010l : iPool ev 2010; Bower et al 
Krauss et al In Section 2, we present our data 

and describe the radio light curve. In Section 3 



we 



present a new distance estimate to V407 Cyg, based 



on optical Na I D absorption lines. Section 4 considers 
and dismisses the possibility that the radio emission from 
V407 Cyg is dominated by expanding thermal ejecta ex- 
pelled in the nova. In Section 5, we hypothesize that 
the radio emission is produced by the Mira wind, simi- 
lar to the radio-emitting ionized envelopes se en around 
other symbiotic systems ([Seaauist et al.|[l"984[ ). but ion- 
ized by the nova outburst. In Section 6, we predict the 
synchrotron signal from V407 Cyg and show that it was 
significant but hidden behind the absorbing screen of the 
Mira wind. Finally, the first radio detection of V407 Cyg 
in quiescence is presented in Section 7, using archival 
VLA data from 1993. We conclude in Section 8. 



2. THE DATA 

2.1. Observational Method 

At the time of the outburst of V407 Cyg, a new cor- 
relator was being commissioned as part of the VLA up- 
grade. We therefore began monitoring V407 Cyg on 2010 
April 3.7, approximately one month after discovery (2010 
March 10), as some of the newly upgraded facility's first 
observations (Krauss et al. 2010). During each epoch, 
we attempted to observed at L (1-2 GHz), C (4-8 GHz), 
K (19-26 GHz), and Q (40-50 GHz) bands, providing 
coverage across the entire VLA frequency range. Un- 
der program AS1039, we obtained data roughly twice 
a week throughout April and May 2010. The spacing 
of our epochs gradually widened, and monitoring was 
continued until April 2012 via programs lOB-233, llA- 
254, and llB-170. Our monitoring of V407 Cyg has now 
bridged more than an entire VLA configuration cycle, 
beginning in its most compact D configuration. During 
the summer (June-August) of 2011, the VLA was in the 
extended A configuration, and V407 Cyg was spatially 
resolved. The details of this imaging will be presented 
in Mioduszewski et al. (2013, in preparation), but fiux 
densities were difficult to derive during this time at some 
frequencies. 

In a typical epoch, we spent two hours total to ob- 
serve all four bands, providing ^10-15 minutes on source 
at each band. At K and Q bands, we obtained im- 
proved pointing solutions on a nearby calibrator, and 
used J2102-I-4702 for gain calibration. At C band, we 
used J2048-I-4310 for gain calibration, and at L band, 
J2038+5119 and J2052-H3635 were used. The absolute 
flux density scale and bandpasses were calibrated using 
3C48. The earlier data (programs AS1039 and lOB- 
233) were obtained with 256 MHz of bandwidth per 
band, split between two independently tuneable base- 
bands. Later observations (programs llA-254 and IIB- 
170) provided 2 GHz of bandwidth, split between two 
1-GHz-wide basebands. 

We also utilized some early data publicly available in 
the archive under program AL733 (PI Lang), obtained 
on 2010 March 27.9 and 2010 March 29.8 at L and C 
bands. This program observed with one-hour-long blocks 
at each receiver band, and used 256 MHz of bandwidth 
(at L band, this bandwidth was split between 1.4 and 1.9 
GHz, while at C band the basebands were set adjacent 
to one another at 5 GHz). 

Data reduction was carried out using standard routines 
in CASA and AIPS. Gain curves and weather tables (us- 
ing the mean of the seasonal average and the day's mea- 
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surements) were applied. Basic editing and calibration 
were carried out on each epoch and receiver band sep- 
arately. Typically, one iteration of phase-only self cali- 
bration was applied, and the two basebands were imaged 
separately. We measured flux densities of V407 Cyg by 
fitting gaussians to the unresolved point-like source, ex- 
cept during the A-configuration (epoch 26), when two 
Gaussians were fit and summed together to represent the 
higher-frequency flux densities. Calibration errors of 5% 
were assumed at L and C bands, 10% at K band, and 20% 
at Q band (Q band was particularly plagued by poor re- 
ceiver performance, a faint gain calibrator, and inclement 
weather) . Measurements are presented in Table [TJ 

Additional r a dio continuum data were pub lis hed by 
iNestoras et all (l20Tol). i Gawronski etldl pOTot ). iPoole'vl 
(|2010D . and IBower et a l. (2010) in Astronomer's Tele- 
grams. Unfortunately, the facilities and frequencies used 
for these measurements are heterogeneous, so they are 
difficult to integrate into our VLA light curve and are 
not used in our analysis (although a spectral index from 
the Nestoras et al. data is listed in Table [2]) . 

2.2. Results 

Our multi-frequency radio light curve for V407 Cyg, 
plotted in Figure [TJ shows a gradual rise and decline of 
radio luminosity over the two years following the nova 
outburst. Due to observational complications, there were 
occasionally small shifts in the central frequencies of our 
bands, as shown in Table [1] therefore, for example, the 
"7.4 GHz" light curve actually includes data centered at 
frequencies between 7.0-7.9 GHz. To minimize depen- 
dence on model assumptions, no spectral index correc- 
tion was applied to our light curve plots to account for 
this small difference. 

The first VLA observations were obtained 17 days after 
outburst at only 5 GHz, and V407 Cyg is already radio 
bright at this time (~7 mJy). Overall, the light curve 
rises and declines, with the higher frequencies peaking 
at larger luminosities and declining first. At 45 GHz, 
V407 Cyg is already fading by the time of our first reli- 
able Q-band observation on Day 33. At frequencies be- 
tween 4.9 and 25 GHz, the source continues to brighten 
for a month or two longer, before turning over and de- 
clining. At 1.4 and 1.8 GHz, the flux density rises until 
Day ~70, and then remains roughly constant at this peak 
flux density for another '-^lOO days before gently declin- 
ing. By Day 300, the flux densities at all frequencies 
converge to ^^3.5 mJy and then remain roughly constant 
for the next year, through the time of writing of this 
paper. 

The entire radio evolution of V407 Cyg is relatively 
gentle and slow, showing gradual rises and broad peaks 
in the 1.4-25 GHz light curves. Typical expectations for 
novae attribute the radio emission to expanding thermal 
ejecta, with luminosity inc reasing in proportion to area 
while it is optically thick (jSeaquist fc Bodd 120081 : Sec- 
tion in — so, oc P. However, the rise of the radio 
light curves in V407 Cyg is signiflcantly slower than this, 
scaling as L„ oc t^-^-^-^. 

Even at early times, the spectral index a (defined 
as S„ cx I'") of V407 Cyg is significantly shallower 
than the a = 2 predicted for optica lly-thick thermal 
emission _fe^ -Seaquis t fc Bodd [2001) . The report of 
INestoras et al.l (2010.) and the early VLA data from pro- 



gram AL733 indicate that the spectrum may be steepest 
at early times, with a > 1.0 (see Table |2|). By the time 
our full VLA monitoring effort began on Day 24, the 
spectral index was measured to be a = 0.7 — 0.8. 

On Day 24, the spectrum could be fit with a single 
power law of a = 0.75 between 1.4 and 49 GHz (Figure 
12) . This spectral characterization remained roughly con- 
stant through Day '^70, with some evidence for gentle 
flattening with time to a « 0.6. After Day 70, the 25 
GHz light curve clearly turns over and the spectra show 
flattening above ~20 GHz. By Day 100, the spectral flat- 
tening has cascaded down to ^6 GHz, and by Day 200 
the entire spectrum is consistent with a flatter power 
law of a = 0.25. At late times (> 350 days after out- 
burst), the spectral index is even flatter, consistent with 
an optically-thin thermal spectrum {a = —0.1; Figure 

ID. 

3. DISTANCE TO V407 CYG 

Our interpretation of the radio light curve de- 
pends on the distance to V407 Cyg, which is poorly 
known. Es timates in the literatur e range from 1 .7 kpc 
(|Kolotilov et al. 200l) to 2.7 kpc (|Munari et al.| [r990). 
all using the period-luminosity relation for Miras. How- 
ever, the period-luminosity relation is n ot constrained for 
the v ery long period of V407 Cyg (e.g., iWhitelock et al.l 
l2008f l. and the relation must be widely extrap olated to 
derive a distance. In addition. iTatarnikova et a l. (2003a]) 
identified lithium in the atmosphere of V407 Cyg. Li- 
rich AGB stars have been noted to lie above the period- 
luminosity relation, displaying anomalously h i gh lu - 
minosities for a given period (|Whitelock et al.ll2003[ ). 
Therefore, distances to V407 Cyg derived from the 
period-luminosity relationship are highly uncertain and 
likely underestimate the true distance. 

To provide additional constraints on the distance us- 
ing interstellar Na I D absor ption lines, we revisi ted the 
echelle spectra published by iMunari et al.l (|2011| ). Two 
velocity components were observed that were constant in 
time as the 2010 outburst evolved (Figure 12]), indicating 
an interstellar foreground origin. They have heliocentric 
velocities of -45.7 ± 1.1 km s"^ and -10.9 ±1.0 km s~^ 
In the local standard of rest (LSR), the corresponding ve- 
locities are —29.9 km s~^ and -1-4.9 km s~^ (for Galactic 
coordinates I = 86.98, b = -0.48). 

We use the Galactic rotation curve of iBrand fc Blitj 
(11993') to translate the LSR velocities of the absorbing 
components to distances of 3.0 kpc and 1.5 kpc, respec- 
tively. Both components are strong and nearly saturated 
(Figure |3]) , indicating a significant amount of absorbing 
material. If the Milky Way's neutral ISM is approxi- 
mated as a slab extending 150 pc from the midplane, our 
line of sight will run through it for ^^3.5 kpc, consistent 
with the existence of absorption from a distance of 3.0 
kpc. This line of sight crosses the Orion Spur spiral arm 
at a distance of ^^1.5 kpc and me ets the large Perse us spi- 
ral arm at a distance of 3 kpc (|Churchwell et all 12009). 
Therefore, the interstellar Na I D absorption implies a 
distance >3.0 kpc for V407 Cyg. 

Several other authors have hypothesized that the more 
blue-shifted absorption component (wlsr = —30 km s""'^) 
is due to the circumstellar envelope of V40 7 Cyg itself 
(ITatarnikova et al.ll2003H iShore et al.ll201lD . If this hy- 
pothesis were accurate, we would only be able to con- 
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Figure 2. Selected radio spectra siiowing tiie evolution of V407 Cyg. The 1-45 GHz power law fits described in Tableware marked as 
dashed red lines. In the bottom-right panel, comparison spectra are shown for three physically significant spectral indices. 

distance would affect our results. 

4. THERMAL NOVA EJECTA AS A POTENTIAL SOURCE 
OF RADIO EMISSION FROM V407 CYG 

Here, we compare the multi-frequency radio light curve 
of V407 Cyg with a model that has historically explained 
the radio emission from most classical novae: expanding, 
thermal cjecta spherically expelled by the nova. We show 
that if this mechanism is the dominant source of radio 
luminosity, the nova ejection in V407 Cyg must have been 
massive (10~^ — 10"'* Mq) — in conflict with significant 
observational evidence that the ejecta were low in mass 
(10-'^ - 10"^ Mq). 

In the past, most measurements of radio emission from 
novae were consistent with thermal bremsstrahlung from 
a ~10* K ionized sphere that is expanding in size, drop- 
ping in density, and gradually tra nsitioning from opti- 
cally thick to optically thin (e.g., iSeaauist fc Palimakal 
[19771: IHiellming et"alT[T979l) . The radio emission has 
been interpreted as the expanding nova ejecta, and 
from basic assumptions about its density profile and 
from simple modeling of its ex pansion, the ej ecta mass 
and energy can be estimated (jSeaauist fc Bo de 2008). 
This scenario has mostly been applied to classical no- 
vae with minimal circumbinary material surrounding 
them, and indeed, the radio emission from the only 
other well-studied nova in a symbiotic system, RS Oph, 



strain the distance of V407 Cyg to > 1.5 kpc. Tatarnikova 
et al. point out that the more blue-shifted Na I D com- 
ponent has a very similar radial velocity as the Ca I 
and Li I absorption features, which presumably arise in 
the Mira's atmosphere. However, we observe that the 
more blue-shifted absorption component did not vary 
in strength relative to the less blue-shifted component 
over the course of the nova outburst, implying that both 
absorption components originate in the interstellar — not 
the circumstellar — medium. In addition, while the nar- 
row absorption component of Ha varies in its central 
wavelength throughout the nova outburst (it shifts blue- 
ward with time), the Na I absorption is remarkably con- 
stant in its central wavelength (Munari et al. 2012, in 
preparation). We therefore conclude that the correspon- 
dence in velocity between the ulsr — ~30 km s~* Na I 
D absorption component and stellar features is a coin- 
cidence, and both Na I D absorption components are 
interstellar. 

Hence, we find that V407 Cyg is located >3.0 kpc away. 
We will take a distance at the lower end of this range, 
-D = 3.0 kpc, throughout this paper, as it is consistent 
with most other works on the 2010 nova, which assume 
a distance of 2.7 kpc. However, we maintain that the 
distance to V407 Cyg is still uncertain. Where appropri- 
ate throughout our work, we note how a change in the 
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Figure 3. Detail of echelle spectra around the Na I D absorp- 
tion lines. We plot eight distinct spectra, observed at a range of 
times after the nova outburst, with the absorption lines superim- 
posed on the time-variable red wing of the He I A5976 line. The 
Na I D absorption in all spectra shows a constant two-component 
structure. 



has proven much more complex (iHiellming et al.l 119861: 
■O'Brien ct al. 2006; Kantharia et al. 2007: R upen et al l 
ToO&.Sokoloski et al.„200S: .Evres et al. 2009.) . Still, here 



we begin with the simplest explanation for radio emission 
in V407 Cyg, and compare our VLA light curve with ex- 
pectations for the expanding nova ejecta alone. 

To first order, most radio light curves of classical novae 
can be fit with a simple model wherein the ejecta are in- 
stantaneously ejected and homologously expanding. The 
model requires a maximum and minimum velocity (wmax 
and Wniin) as input for this "Hubble flow" , and the ejecta 
are assumed to be distributed in a thick expanding shell 



with a p cx 



radial density profile. In this simple 



model, the ejecta are also assumed to be isothermal, with 
an electron temperature ~ 10^ K. Additional complica- 
tions can certainly be added to this model, like clumping, 
temperature increases, and non-spherical geometry (e.g., 
EjcUming 1996; Hcywood 2004; Rov ct al. 2012). How- 
ever, in most cases, the simplest model has been sufficient 
to approximate observed light curves. 

In Figure m we superimpose such models on our VLA 
light curve for V407 Cyg to test if the ejecta might 
account for the majority of the radio emission. We 
assume a distance of 3.0 kpc, Vmax = 3,200 km s~^ 
(|Munari et al.ll2011l: iShore et al. ] |2011[ ), and a thick sheU 
with Vinin/^'max = 0.1. Although no Hubble- fiow model 
provides a satisfactory fit to the data at all frequencies, 
the top panel shows a relatively good fit to the higher 
frequencies and earlier times. This model requires a 
slightly higher electron temperature than typically as- 
sumed, Toj ~ 4 X 10'' K (alth ough consistent w ith ob- 
servations of some novae; e.g., lEvres et al.lll996[ ). to fit 
the high fiux densities and relatively early peak at 41/45 
GHz. It also requires a massive ejection (Moj « 1 x 10~* 
M0) to achieve the high measured radio fiux densities. 



Such a large ejecta mass is typical of classical novae, 
but contradicts optical and X-ray observations which im- 
ply that the nova in V407 Cyg was a relatively low-mass 
ejection, more typical of recurrent novae than classical 
novae. The optical light curve decl ined by three magni- 
tudes in the V-band over 24 days (jMunari et al.l[2011.' ). 
According to nova models, s uch a rapid decline should 
only o ccur if Me\ < 1 "^ M f:, (jYaron et al.ll2005D . In ad- 
dition, iMikolaiews^ (|2010[) points out that V407 Cyg's 
light curve matches the light curve of RS Oph's 2006 
nova event remarkably well; estimates of the ejecta 
mass for RS Oph r ange across ~ 10~^ — 10~^ Mp 



(lO'Brien et al.l 119921: [Sokoloski et al 



120091: lEvres et al.ll20Q9l: iVavtet et a! 



20061 iDrake et al l 
20m. Additional 



support for a low ejecta mass is lent by a soft X- 
ray component, observed in V407 Cyg between Days 

15-50 and identified as thermal emissi on from resid- 
ual nuclear burning on t he WD surface (jSchwarz et al.l 
120111: iNelson et aLri2012D . The speed with which the 
ejecta become optically thin is also confirmed by optical 
spectroscopy, where high-ionization lines (e.g., [Fe XI]) 
first appear around Day 7 and subsequently brighten 
(Munari et al. 2010,). The rapid "turn-on" time of soft 
X-ray emission and coronal lines imply a low ejecta mass 
of ~ 10~^ — 10"'' M(7), assuming an ejecta velocity of 
3,200 km s-i (ISchwarz et al.ll20Tl . 

A small ejecta mass of Mcj = 1 x 10~^ Mq yields the 
model shown in the bottom panel of FigureHl This model 
underpredicts the radio fiux densities by two orders of 
magnitude. The model light curves also evolve extremely 
quickly, with the predicted 5 GHz curve peaking ~ 20 
days after the start of outburst rather than at Day ^ 100 
as observed. 

Even if the ejection in V407 Cyg was massive, seri- 
ous problems remain in fitting the radio spectral energy 
distribution with the model pictured in the top panel of 
Figure 21. The figure illustrates this point; the model 
matches the 45 GHz light curve rather well up until 
Day ~160, but offsets between data and models grow 
toward lower frequencies. At 1.8 GHz, there is roughly 
an order of magnitude discrepancy between the observed 
and model light curves. This discrepancy is caused by 
the model prediction of an optically-thick thermal spec- 
trum at early times (a = 2), which is much steeper 
than the a = 0.6 - 1.0 that we observe (Table E]). The 
model also predicts that, once the thermal emission be- 
comes optically thin, it will rapidly fade. In contrast, 
we have observed a late-time plateau in the light curve 
of V407 Cyg with a duration of >300 days (continuing 
through present) during the time when the spectral index 
implies optically-thin emission. 

The simple thermal model assumes that the ejecta are 
not clumpy, but smooth in density with a volume fill- 
ing factor of / = 1. However, observations of novae 
reveal that the eje cta can be quite clumpy, with / = 
10"^ - 10"^ (e.g., iS aizar fc Ferland 1^9^ MasoneLaLl 
l2005HEderoclite et a l. 2006; Shara ct al. 2012). In a sim- 
ple model, for isothermal ejecta with a single filling fac- 
tor describing all times and radii, dumpiness simply in- 
creases the radio luminosity by a constant factor over the 
uniform-density case, wi th no changes to th e time evolu- 
tion or radio spectru m ([Abbott et al.lll98"lL iNugis et"all 
119981: iHevwoodI 120041 ). This basic case is likely an over- 
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Figure 4. Models for the multi-frequency radio emission produced by thermal expanding nova ejecta are superimposed as dotted lines on 
our VLA light curves for V407 Cyg. The models assume a "Hubble flow" distribution of ejecta, a maximum velocity of 3,200 km s~^, a 
distance of 3.0 kpc, and a ratio of fmax/''min = 0.1. The top panel shows a model for a massive ejection that fits the data relatively well 
(Moj = 1 X 10^* M0, Tq = 4 X 10^ K), while the bottom panel shows a model consistent with multi-wavelength observations of V407 Cyg 
and expectations for recurrent novae (Mej = 1 X 10~® M©, Tgj = 1 X 10* K). Note that the vertical scales of the top and bottom plots 
differ. 



simplification, and more complex modelling of the ra- 
dio emission from multi-temperature and time-variable 
clumped ejecta is needed. However, it is difficult to imag- 
ine that the addition of such complications could slow the 
evolution of the low-mass model light curve by the factor 
of ~4 needed to match the data. Therefore, while clump- 
ing might serve to increase the radio luminosity from our 
Mej = 1 X 10~^ M0 model to bring it into better align- 
ment with observations, clumping can not account for 
the dramatic discrepancies in light curve timescales and 
spectral indices. 

Of course, the nova ejecta in V407 Cyg are decelerat- 
ing with time as they interact with the wind from the 
Mira giant companion (clearly seen in opt ical line pro- 
files; iMunari et all 120 lit iShore et al.ll20lH) . In the pro- 
cess of decelerating, the outer ejecta will be swept up into 
a relatively thin and dense shell; the radio emission will 
therefore remain optically thick longer than in the con- 
stant expansion model, and then transition more quickly 
to an optically thin state. Such an evolution would only 



exacerbate the disagreement we already see between the 
data and models of expanding nova ejecta — -namely, an 
observed spectral index at early times which is much flat- 
ter than a = 2 and a more-gradual-than expected rise 
and decline in the observed light curve. This consider- 
ation serves as a reminder that significant circumbinary 
material surrounds the nova, which may itself contribute 
to the radio emission (see the next Section). 

The required ejecta mass shrinks, but remains quite 
high (^4xl0~^ M0) if the distance is small, ~1.5 kpc. 
Independent of distance and clumping, a model with 
^ 10~^ Mq peaks far too quickly and declines 
too rapidly, compared with our observed light curve for 
V407 Cyg. In view of the poor fits to expanding thermal 
ejecta and the well-established presence of circumbinary 
material around V407 Cyg, we now proceed to consider 
other physical explanations for the radio emission from 
V407 Cyg's 2010 nova explosion. 

5. RADIO EMISSION FROM AN IONIZED MIRA WIND, 
SHOCKED FROM WITHIN 
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Figure 5. The set-up for our model of V407 Cyg. The size of the blue and red circles, denoting the WD and Mira giant respectively, are 
larger than to-scale for clarity. A wind emanates from the Mira giant, depicted as an orange radial gradient. Along a given polar angle 
9, the radial density profile of the circumstellar material is approximated as independent of azimuthal angle 0. The radial density profiles 
depend on the separation between the WD and Mira giant (a) and the wind mass-loss parameters and t)„ . Two such density profiles 
are plotted, for 9 = ±20° and 9 = ±170°, assuming = 2 X 10~^ Mq yr'^, Vw = 10 km s'^, and a = 17 AU. We calculate r = 1 
surfaces and determine their radii from the WD at each of an array of impact parameters (z), assuming that the system is being viewed in 
the orbital plane, with the Mira giant between the viewer and the WD (from the right, in this figure). 



We now consider the contribution of the Mira wind 
to the radio emission from V407 Cyg. The ionized en- 
velopes surrounding symbiotic stars emit thermal radio 
emission; their luminosities scale with the wind density 
and the degree of ionization, but span the range of radio 
luminosities observed fro m V407 Cyg after the 2010 nova 
(jSeaquist k Taylori[l990l ). We show that the Mira wind 
in V407 Cyg can produce the bulk of the radio luminosity, 
as it becomes ionized by the hot WD. The radio lumi- 
nosity grows while the wind is becoming more ionized, 
and declines as this ionized wind gradually recombines 
and becomes swept up by the nova shock. 

5.1. A Quantitative Model 

Here, we expand on the simple model introduced in 
iNelson et all (|2012D , which was originally intended to de- 
scribe the evolution of the nova blastwave in the asym- 



metric circumbinary medium surrounding the WD in 
V407 Cyg. We solve for the ionization of the wind as 
a function of radius, polar angle from the WD, and time. 
We compare the location of the ionization front with the 
location of the nova shock front (also a function of polar 
angle from the WD and time), to estimate the portion of 
the wind that has been shock-heated and no longer emits 
in the radio band. Then, assuming a viewing direction, 
we can calculate the radio flux density and model the 
radio ligh t curve. A simila r mechanism has been sug- 
gested by iKennv &: Tavlorl ([2005) to explain the radio 
flux densities and images of symbiotic stars in outburst; 
they model these systems as an ionized, fast, hot wind 
emanating from the WD, colliding with a slower, cooler, 
partially-neutral wind from the red giant. 
V407 Cyg is an inherently asymmetric system; photons 



emanating from the WD in the direction of the giant en- 
counter significantly denser material that those traveling 
away from the red giant. Our model accounts for this 
asymmetry by splitting the environment of the binary 
into slices of polar angle (9) centered on the WD posi- 
tion (Figure [5|). Assuming spherical symmetry around 
the Mira giant, a given slice can be described with a sin- 
gle radial profile of circumstellar density, independent of 
azimuthal angle (</>) and time. The density profile in a 
given polar-angle slice depends only on the binary sep- 
aration a, the mass- loss rate of the wind M^, and the 
wind velocity Vuj- Relative to the Mira giant, the density 
of the spherically symmetric circumbinary material is: 



47r /im^fWu 



' RG 



(1) 



where ncsM is the number density of particles in the 
wind at a radius trq from the Mira giant; fj, is the 
mean molecular weight and mn is the hydrogen mass. 
Throughout this paper, we quote assuming — 10 
km s~^, as e stimated for V407 Cyg from optical P Cygni 
line profiles (jAbdo et al.ll2010D . Relative to the WD, the 
circumbinary medium is asymmetric as a function of 6: 



ncsMirwo) 



[{a - rwD cos6')^ -f- (rwD sinfl)^] 



(2) 



where 6* = is the direction pointing from the WD to- 
wards the Mira giant and rwD is the distance from the 
WD. Radial density profiles for two different values of 
9 (one toward the Mira giant and the other away) are 
plotted in Figure [5] Regions of the Mira wind away from 
the giant will become ionized more quickly than toward 
the giant (illustrated by comparison of our model's den- 
sity profiles for the 9 = ±170° and 9 = ±20° directions, 
shown in Figure [5]). Also, the nova blastwave will expe- 
rience less deceleration in the = ±170° direction than 
in the 9 = ±20° direction. 

5.1.1. Photoionization Model 

We assume that, at the time of nova outburst, the 
Mira wind is completely neutral. We then input a time- 
variable flux of ionizing photons {Nu(t), in units of pho- 
tons s~^), assumed to be emanating from the WD as a 
result of nuclear burning on its surface. Our selections for 
this important input are described more below, in Sec- 
tion 15.1.41 The fraction of the ionizing flux that enters a 
given polar-angle slice centered at 9 is; 



. 1 
/vol = 2*^°^ 



A9\ 1 / A9\ 

-cos [9 + — ]. (3) 



The width of the polar-angle slice is A9, and set to 1 deg 
in our model (for a total of 180 slices). 

At a given time step, a given polar angle slice is split 
into radial bins of width Ar; starting at small rwo a-nd 
traveling outward, we calculate how much ionizing flux 
reaches a given radial bin and the ionization fraction 
x{r,9,t) in that bin,. The ionizing flux at the WD pho- 
tosphere is Su{r = 0, 9, t) = /vol Nu{t). The ionizing flux 
reaching radial bin (r, 9, t) is Su{r — 1, 9, t), while the ion- 
izing flux escaping bin (r, 9, t), after ionizing the neutral 



atoms in this bin, is Su{r, 9, t). The relationship between 
incident and emerging radiation is: 

Su{r,9,t)^Su{r-l,9,t) 

"4^ [1 - x{r, 9, t - 1)] nc.„,(r, 9) r^^ Ar /voi(0) 



At 



(4) 



- ['iTTx{r,9,t - l)'^ ncsm{r,9)'^ ai^2) r^jD Ar f^oii 



All parameters are in cgs units. The second term on the 
right hand side represents the number of neutral atoms 
in this radial bin available for absorbing the ionizing flux. 
The third term is the rate of recombinations in this bin — 
newly neutral atoms that will also absorb Su{r — l,9,t). 
Here, a(2) is the recombination coefficient of hydrogen, 
excluding recombinations to the n = 1 level. At a tem- 
perat ure of 10^ K, a(2) = 2.76 x lO'^^ g-i (jSpitzeii 
Il978f ). This value is conservative, as recombination plays 
a relatively weak role in the evolution of radio emission 
from V407 Cyg; hotter assumed temperatures would give 
lower values of q;(2) a-^d even weaker effects from recom- 
bination. 

We assume that the ionization front is sharp; all atoms 
in a bin must become ionized {x{r,9,t) = 1) before the 
residual ionizing flux is allowed to escape to the next ra- 
dial bin. In support of this assumption, we only include 
photons with energies between the Lyman limit and 5 
times the Lyman limit (13.6-66.5 eV) in the count of 
ionizing photons. The hardest photons are excluded due 
to their small cross sections for ionizing hydrogen; they 
presumably serve to ionize metals, although we include 
no treatment of metals in our simple model. If the ioniz- 
ing flux reaching a bin is not sufficient to ionize all of the 
atoms in this bin (in other words, if Su{r,9,t) in Equa- 
tion 2] is negative), the fractional ionization of this bin is 
calculated as: 



x{r,9,t) = x{r,9,t- 1) + 



Suir~l,9,t) At 



Ait ricsmir, 9) r"^^ Ar /vol (6*) 
47r x(r, 9, t - 1)^ ncsm(r, 9 f a(2) r^^ Ar ^01(6*) At 



47r ncsm{r, 9) Ar /vol (6*) 

(5) 

This formulation properly accounts for those time steps 
and bins which are not completely neutral {x > 0) and 
are not reached by any ionizing flux; these will slowly 
recombine in the absence of an ionizing radiation field. 

5.1.2. Nova Shock Front Model 

Because 10^ K ionized gas has very different radio- 
emitting properties than 10^ K gas, we must also track 
the location of the nova shock as a function of time and 
9. Increasing the electron temperature by three orders of 
magnitude decreases the radio optical depth by four or- 
ders of magnitude (Equation[B]). The nova blastwave will 
sweep up circumbinary material and heat it from radio- 
emitting temperatures to X-ray emitting temperatures. 

Our model of the evolving nova blastwave is as de- 
scribed in Ne lson et al.l (|2012t ). We input M^, v^, a, 
Afcj, and initial blast wave velocity. The model assumes 
that the ejecta are expanding in a thin shell at a single 
velocity initially. Each polar-angle ring then has an in- 
dependent evolution of radius and velocity that depends 
on the density profile it encounters. It decelerates as it 
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sweeps up material, conserving kinetic energy in a Sedov- 
like evolution. 



5.1.3. Calculation of Radio Flux Densities 

The luminosity of thermal radio emission depends crit- 
ically on the emission measure (EM), here defined as 
J nl^dl, where is the electron density of the absorbing 
material and dl is the path length through this material. 
The optical depth Ti, at an observing frequency v along 
a given line of sight is: 



and sum them up. The EM for each bin is calculated as: 



/7;\ 1-35 , y .-2.1 / EM 

= ^-^'^ ^ j (gH^) 

, (6) 

(|Seaauist fc Bodd 120081) . The flux density produced 
by this line of sight is: 



(^)-— (^) i^y (I 



(l-e--) 



GHz 



(7) 



Here, A is the projected area of this line of sight, D is 
the distance to the emitting source, and is its electron 
temperature. 

In our model of V407 Cyg described above, we can 
calculate S^, if we assume a viewing angle; optical spec- 
troscopy observationally constrains this an gle fo r the 
2010 outburst of V407 Cyg. iMunari et all (|20lH) and 
I Shore et all (|2012^ point out that narrow absorption lines 
of hydrogen and metals are present in the spectrum of 
V407 Cyg following the nova outburst, implying the pres- 
ence of neutral circumstellar material. Although these 
lines fade with time, they never completely disappear, 
implying that the wind never becomes completely ion- 
ized. Therefore, the nova was probably observed while 
the giant was in fr ont of the WD along our line of sight 
([Shore et al.1 120121) . Our best guess at the o rbital pe- 
riod of V407 Cyg (43 years; IMunari et allll990[ ) is signif- 
icantly longer than the timescale of our observations, so 
we assume this viewing angle is constant across our radio 
light curve. We show a schematic picture of the binary 
in Figure [5l which we assume is being viewed from the 
right-hand side of the figure. 

While our model is not currently equipped to provide 
light curves for an arbitrary viewing angle, the radio light 
curves presented here should be relatively unaffected by 
small changes in viewing angle, as the radio luminos- 
ity is predominately dependent on the projected angular 
area of warm ionized wind, which will weakly depend on 
viewing angle. This insensitivity to viewin g angle is illus- 
trated for symbiotic stars in Figure 6 of ISeaauist et al.l 
([1984), where, for winds that are significantly ionized, 
the measured flux densities and spectral indices change 
negligibly with a 90° change in viewing angle. 

We assume an electron temperature and construct an 
array of impact parameters (values of separated by 
Az), along the axis perpendicular to our line of sight 
(Figure [5]). For a given impact parameter, we find all 
bins of (r, 6*) that are intercepted by this line of sight. 



^ , ^ N x(r,9,t)'^ ncBm(r,9Y 
EM{r,0,t) = ^ \\L X 



3.086 X 1018 

z 



itm{e 



(8) 



cm ^ pc. 



However, a bin is only considered as contributing to 
the radio-emitting EM if it has not been shocked by 
the nova blastwave. If a bin has been swept over by 
the nova shock, c a lculat ed as described in Section [5.1.21 
and iNelson et al.1 ()2012[ ). it will typically be heated to 
10^ — 10^ K, and will therefore be optically thin and 
faint at radio wavelengths; we set EM(r, 6*,^) = in 
this case. We calculate A by assuming that the emit- 
ting area at a given impact parameter is ring-shaped: 
A = iT{{z + ^f-{z-^f). 

To understand the evolution of the optically thick ra- 
dio emission, we find it instructive to calculate surfaces 
where the optical depth reaches a value of unity, as a 
function of time. Radio observations at a given frequency 
V will be insensitive to conditions within the = 1 sur- 
face. The size of a = 1 surface depends on the path 
length through and density of the ionized gas. To reach 
— 1, an EM of ionized material must be obtained 
equivalent to: 



EM^ 



pc 



12.1 



1.35 

(cife) 



2.1 



(9) 



To calculate the shape and size of the — \ surface 
at a given time and frequency, we again step through 
values of z. At a given impact parameter, we proceed 
from the bin closest to us along our line of sight, located 
far out in the Mira wind, inward toward the binary and 
potentially beyond it, summing up the EM in bins until 
EMt- is reached. The bin where EM,- is obtained marks 
a location in (r, 9) of the r = 1 surface. This process is 
iterated for all impact parameters to obtain an array of 
such locations and map the topology of the r = 1 surface 
(plotted as black lines in Figure E]) . 

In Figure [71 we plot the evolution of the radius of the 
r = 1 surface at three different frequencies, measured 
perpendicular to our line of sight. The radius of the 
r = 1 surface is defined as the maximum impact param- 
eter where r = 1 is reached. In some models, the r = 1 
surface does not always look like a filled circle, but dis- 
plays a ring-like geometry at later times, because EM^ 
is never reached at small values of z (e.g., Day 800 in 
Figure [6]). Ring-shaped t — 1 surfaces are produced by 
the nova blast eating away at the thermal wind from the 
inside. They are depicted in Figure [7| as having both a 
minimum and maximum radius (e.g., the 1.4 GHz curve 
in Figure [71 between Days 500 and 1000). 

5.1.4. Assumed Physical Parameters 

In exploring the parameter space of V407 Cyg, we hold 
most parameters fixed and focus on varying the wind 
mass loss rate, AI^,. The wind velocity is assumed to be 
10 km s^i in our discussion of M^. The binary separa- 
tion is held constant at 17 AU. In calculating the evo- 
lution of the nova blastwave, we assume an initial shock 
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Figure 6. Models of V407 Cyg, with the binary located at the origin (small red and blue dots; not to scale) and the nova expanding 
into a Mira wind of mass-loss rate = 2 X 10~^ yr"-"^, wind velocity Vn, = 10 km s~^, and wind temperature Te,w = 1 X 10* K 
(orange gradient), as pictured in Figure [5] We show the environment of the binary at four different epochs: 1, 25, 200, and 800 days after 
outburst. In each epoch, the portion of the wind that is ionized is filled with white and outlined in blue. The extent of the nova shock 
front is outlined in green (too small to be visible here on Days 1 and 25). Surfaces where optical depth r = 1 is reached are shown as black 
lines, and are depicted at three different frequencies: 1.4 GHz (dashed), 5.0 GHz (solid), and 25.0 GHz (dotted). These surfaces grow as 
the wind becomes more ionized, but are destroyed when the nova shock sweeps over them. 

V407 Cyg inform our assuraptions of how the blackbody 
teraperature and luminosity vary with time. We plot the 
assumed evolution of the hot WD in the bottom three 
panels of Figure [T] 

At the onset of the thermonuclear runaway, a brief 
"flash" of ionizing radiation is released, before the ejecta 
has time to expand and reprocess the hot WD spectrum 
to a puffed-up ~10,000 K blackbody. We assume this 



veloc ity of 3,200 km s"^ (jMunari et al.ll201lHShore et al.l 
1201 ID and an ejecta mass of 10~^ Mq. 

We assume a wind electron temperature of T^^w = 
1 X 10** K, but the calculated radio emission is relatively 
insensitive to Te^w In Figure 13 we plot the same model, 
one with Tg = 1 x 10^ K (solid lines) and the other 
with Te^tu = 3 X 10^ K (dashed lines). At most points in 
the light curve, the models differ by less than a factor of 
~2. 

An essential input for our models is the ionizing flux 
that emanates from the nuclear-burning WD, and how 
this flux varies with time. We assume the WD has a 
blackbody spectrum; theory and X-ray observations of 



flash lasts 1 hour, has a lumin osity of 1.2 x 10'^^ erg 
s~^, and a temperature of 10^ K ()Prialniklll986( ). Subse- 
quently, the ionizing flux dramatically diminishes, until 
the ejecta have expanded sufficiently to become optically 
thin to the hot WD once more. 
There is evidence for significant ionizing flux in the 
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Figure 7. Panel (a): For the same model parameters as depicted 
in Figure [G] we plot the radius of the t = 1 surfaces, perpendicu- 
lar to our line of sight, as a function of time for 1.4 GHz (gold), 
5.0 GHz (red), and 25.0 GHz (blue). The maximum radius of the 
nova blastwave is also shown as a solid black line. Panel (b): The 
assumed evolution of the WD blackbody temperature. Panel (c): 
The assumed evolution of the WD blackbody luminosity. Panel 
(d): The ionizing photon flux produced by the assumed WD tem- 
perature and luminosity evolution. 

Swift X-ray light curve from a s oft X-ray compo- 
nent present between Days ^15-50 (iShore et al.l 120111 : 
ISchwarz et alllMlt INelson et al.ll2012D . The soft X-ray 
component is likely attributable to the hot WD, which 
would then be the primary source of ionizing flux in the 
V407 Cyg 2010 nova. Note that the hard X-ray compo- 
nent (attributed to forward-shocked material in the nova 
blastwave) is observed to be brighter than the soft X-ray 
component in the Swift lightcurve, but this is only an 
observational effect produced by the severe absorption 
of the soft X-rays by circumstellar material. Because of 
the strong absorption, it is difficult to precisely deter- 
mine the luminosity of the hot WD from the Swift data 
(jSchwarz et all 120111: iNelson et al.ll2012D . We therefore 
assume a WD plateau luminosity of 1 x 10^^ erg s~^, 
as consistently implied by mode ls of nuclear burning o n 
WDs during nova events (e.g., iSala fc Hernanzl [20051 ). 
We do use the Swift measurements to constrain the rising 
temperature of the WD up until a peak temperature on 
Day 39. The rise is fit roughly with Two = 2 x 10* t^-^'^ 
K, where t is the time since outburst in units of days. 
Note that this fit is extrapolated back to the time of on- 
set of outburst, and implies that the ionizing flux from 
the WD is peaking just a few days after outburst (Fig- 
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Figure 8. Light curves for V407 Cyg at three frequencies: 1.4 
GHz (gold), 5.0 GHz (red), and 25.0 Ghz (blue). The diamond- 
shaped points connected by dotted lines are our VLA data. Models 
assuming A/,„ = 2 X 10~® Mq yr"^ and Te,™ = 1 X 10"' K (as in 
Figure [6| are shown as solid lines. Dashed lines arc models for 
M„ = 2 X 10-8 Mq yr-i and T^,™ = 3 X 10"* K. 

ure [7]). This early-time evolution is highly uncertain 
and unconstrained by data, and yet unfortunately it is 
the primary determinant of the early time radio light 
curve. After the temperature peak on Day 39, we as- 
sume the luminosity of the WD declines from its plateau 
value as Lwd oc t"^ -'^*, with a roughly constant radius 
so TwD oc (|Prialniklll986l) . 

5.2. Interpreting the Light Curve with Our Model 

The ionization state of the Mira wind varies in the 
aftermath of the nova explosion, as evidenced by op- 
tical spectroscopy. Narrow Ha emission lines are seen 
at early times (just tw o days after outbur st), consistent 
with an ionized wind ()Munari et all 1201 It ). In the first 
months after explosion, narrow P Cygni profiles were ob- 
served for low-ionization metal lines (Shore et al. 2012). 
These P Cygni metal lines disappeared by Day 70, while 
the Balmer line a bsorption weakened over this period 
(jShore et al.ll2012[ ). Shore et al. attribute this to increas- 
ing ionization in the Mira wind over the months follow- 
ing outburst. This increase in ionization of the wind 
is not surprising, given observations of a hot WD be- 
tween Days ~ 15 — 50 (Section l5.1.4|) . It is worth noting 
that some Ha absorption was sustained at all times, im- 
plying that the win d never becomes completely ionized 
ijMunari et"aI|[20TTl ). 

In our VLA data, the ionized wind of the Mira giant 
radiates as free-free emission, just as detected for other 
symb iotic stars in quiescence (e.g., iSeaquist fc Tavloil 
ll99Clf ). The spectral index at early times (a = 0.6 — 1.0) 
is consistent with the theor etical expectation for a spher- 
ical p oc r" wind (g = q.6:l01non| [l975HPanagia fc Fellil 
119751 : iWright fc Barlow! 119751 ) arid spectral index mea- 
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surements for symbiotic stars (a ranging from — 0-1 to 
1.4. with a median of 0.8: lSeaauist fc Tavloiill990D . How- 
ever, V407 Cyg is more complex than the typical sym- 
biotic star, because the radio flux densities are changing 
rapidly with time. Our model shows that the 1.4-25 GHz 
hght curves rise at early times because the wind becomes 
more ionized (Figures E| and |H1) ■ In other wo rds, the X 
parameter, as defined in lSeaauist et al.l ()1984D to param- 
eterize the fractional ionization of the wind, is increasing 
with time. Seaquist et al. show that as X increases, the 
spectral index of the partially optically-thick emission 
flattens — consistent with the decrease in a from 1.0 to 
0.7 between Days 12 and 39. 

As X increases and the "cone" of neutral wind de- 
creases in extent, the ionized emission measure at the 
outer extents of the wind increases. The radio flux den- 
sity thereby brightens and the t^, = I surfaces grow in 
diameter (e.g., compare Days 1 and 25 in Figure IH]). We 
stress that our radio light curve can not probe condi- 
tions internal to the Ti, sa 1 surface. These surfaces 
are at radii of > 10^^ cm, significantly larger than the 
nova blastwave radius (black solid line in the top panel of 
Figure ^ . Complex processes may be occurring within 
the — 1 surfaces, like thermal emission from expand- 
ing nova ejecta (Section 2]) or synchrotron emission from 
blastwave interaction (Section 15]) . However, this emis- 
sion would only become apparent when either the wind 
becomes optically thin or the additional emitting com- 
ponents expand to radii larger than the w 1 surface. 

As can be seen in Figure [U the model agrees with our 
observed radio light curves as well as can be expected, 
given the significant uncertainties in the ionizing source. 
The model plotted in Figure[8]assumes a Mira wind mass- 
loss rate of = 2 x 10~^ Mq yr~^, a distance of 3 kpc, 
and Te^w = 10'' K. It broadly reproduces the observed 
features of the radio light curve. The peak flux densities 
are roughly matched (to within a factor of ~2; compare 
this with the order-of-magnitude discrepancy at lower 
frequencies in Figure |4]), as are the slow decline times 
and the earlier turn-over at higher frequency. 

However, the model light curves peak too early, on 
Day ~15, relative to the observed maxima of Day ~ 
30 — 100. Revising the input ionizing radiation field 
and its evolution — a tweak completely justified given the 
paucity of observations constraining the hot WD at early 
times — would ease this discrepancy. We have no con- 
straints on the WD's luminosity or temperature before 
Day 15 or after Day 50, while its spectrum is peaking in 
the UV rather than the X-ray — even though the WD is 
likely an important ionizing source at these times. We 
assume that the increase in temperature and luminosity 
of the WD is gradual; extrapolating from the observed 
temperature and luminosity points implies that the ion- 
izing flux is actually peaking very early, around Day 4 
(bottom panel of Figure [?])■ If, instead, the tempera- 
ture and luminosity of the WD incr eased suddenly, as 
observed in RS Oph bv lOsborne et a l. (2011), the pre- 
dicted rise of the radio light curve would change signif- 
icantly. More detailed modeling of the luminosity and 
temperature evolution of the WD is outside the scope of 
this paper. 

If the mass-loss rate is significantly lower than — 
2 X 10"^ Mq yr~i (solid lines in Figure ^ , the fiux den- 
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Figure 9. Light curves for V407 Cyg; similar to Figure \S\ but 
showing models with different M^. Diamond-shaped points con- 
nected by dotted lines are our VLA data. Solid lines show a model 
with Mw = 4 X 10~^ Mq yr"-*^. Dashed lines are models for 
A/u, = 1 X 10~^ M0 yr~l; the 1.4 GHz model is so faint that 
it falls outside the plot limits. Both sets of model lines assume 

Te,w = 1 X 10" K. 

sities never become large enough to match our observa- 
tions, even if the entire wind is ionized. We can under- 
stand this by applying the standard equation for deriving 
mass- loss rates of fully- ionized stellar winds from thermal 
radio flux densities (Wri ght fc Barlow ,197 5): 



M„ = 4.5 X 10^ 
D 



10 km s" 



3/4 



1 mjy / 



3.0 kpc 



3/2 



IGHz 



-1/2 



(10) 



The Gaunt factor gff is a weak function of the wind 
electron temperature: 



= 9.77 -0.17 + 0.130 log 



Te 



3/2' 



(11) 



We estimate Si, from the broad maxima of our 1.4- 
25 GHz light curves, and find A/„ w (6.6 - 7.7) x 
10^*" Mq yr~^ for a distance of 3.0 kpc and a wind 
electron temperature of Te,™ = 10^ K. This estimate 
is only very weakly dependent on the wind tempera- 
ture. For T^,^ = 30,000 K, M„ decreases slightly 
to (6.0 - 7.1) X 10"'^ Mq yr-^ These are essentially 
lower limits on AIw Of course, if V407 Cyg is at a 
smaller distance, the implied Mm will drop accordingly 
{M^ sa 2.5 X lO"'^ Mq yr"! at 1.5 kpc). 

If the mass-loss rate is significantly higher than = 
2 X 10^^ Mq yr~^, as shown by the dashed-line model in 
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Figure IHl the hot WD will only be able to ionize a rela- 
tively small volume. As a result, like in the low case, 
we predict lower-than-observed flux densities and small 
r = 1 surfaces (although, in this case the flux density 
will depend rather sensitively on the assumed luminosity 
and evolution of the ionizing radiation held). 

Not only do small r = 1 surfaces produce too little 
radio luminosity, they also lead to radio light curves that 
decline too quickly, because the nova shock front sweeps 
over them at earlier times. For example, in the high 
Mui case plotted in Figure [HI the model predicts that the 
radio emission at 25 GHz should be optically thin around 
the time the observed light curve is peaking. 

Radio flux densities eventually decline for two reasons: 
recombination and heating by the nova shock. A good 
example of decline by recombination can be seen in the 
1.4 GHz curve in Figure [T] the maximum radius declines 
from 11 X 10^^ cm to ~ 4 x 10^^ because of recombi- 
nation in the outer parts of the wind. At the same time, 
between Days 600 and 1000, we see that the minimum 
radius of the 1.4 GHz r = 1 surface increases from cm 
to ~ 4 X 10^^ cm, until the minimum and maximum radii 
meet and the system becomes completely optically thin 
to 1.4 GHz radiation. 

In our observed light curves, all frequencies stop 
rising after Day ^ 70 — consistent with the findings of 
iShore et al.l (|2012f ) that the wind becomes increasingly 
ionized up through Day ~70. Soon thereafter, the ob- 
served light curves turn over and begin to decline for 
the higher frequencies (>4 GHz), but the 1.4/1.8 GHz 
light curves are constant for another '^100 days before 
declining and becoming optically thin. Correspondingly, 
our model 1.4 GHz curve is quite flat-topped until Day 
'^SOO, when both recombination and the nova blast begin 
to significantly diminish the radio emission measure. 

At late times, the radio flux densities observed at all 
frequencies converge to a steady (or perhaps slowly de- 
clining) '^3.5 mjy (starting around Day 350, continu- 
ing through the time of submission). The radio spec- 
trum is consistent with optically-thin thermal emission 
{a w —0.1; Table[2]). Such a leveling-off of the radio light 
curves is, to date, consistent with our rough model for 
the recombining Mira wind (Figure [8]). Occasional moni- 
toring of V407 Cyg in the upcoming years will show if the 
radio luminosity continues to gradually fade as expected. 



6. WHY NO RADIO SYNCHROTRON EMISSION 
OBSERVED IN V407 CYG? 

From the gamma-ray detection of V407 Cyg, we know 
that particles w ere accelerated to relativistic speeds by 
the nova shock (jAbdo et al.ll201^ . We might therefore 
expect accompanying radio synchrotron emission. How- 
ever, at no time are our measured spectral indices (Ta- 
ble indicative of optically-thin synchrotron emission, 
for which we would expect a = —0.5 to —1.0. Here we 
model the synchrotron emission from V407 Cyg to show 
that bright synchrotron emission was likely produced in 
V407 Cyg, but it was hidden behind the absorbing screen 
of the Mira wind. 

We assume that as the nova shock wave interacts with 
surrounding circumstellar material, it converts a con- 
stant fraction (e) of the post-shock energy density into 



relativistic electrons, such that: 



No E P dE ^ €^1 ruH ncsM v"^ 



(12) 



0.51 MgV 



where all quantities are in cgs units. The blastwave ve- 
locity is Vs and the Mira wind density at the location 
of the blastwave is ncsM (as described in Equation [T]) . 
Here, p is the index of the electron energy spectrum, as- 
sumed to have a power-law form. Abdo et al. (2010) es- 
timate p with large uncertainties from the gamma-ray 
spectrum: p — 1.75lo;g. We can also directly mea- 
sure p if optically-thin synchrotron emission is detected, 
as it is in a few other novae. The radio spect ral in- 
dex im plies p = 2.1 — 2. 5 in RS Oph (Rupen ~et al.l 
'2008: Sokolosk i et all [20081) and p = 2.4 in GK Persei 
(Seaquist et al.lll989f ): we assume p = 2.3 for V407 Cyg. 
The minimum energy in the integral, 0.51 MeV, corre- 
sponds to the rest mass energy of the electron and ap- 
proximates the energy at which electrons become rela- 
tivistic. The efficiency of parti cle acceleration, e, is as- 
sumed to be constant with time. lAbdo et al.l (|2010f ) show 
that, if the gamma-ray emission was solely produced by 
leptonic processes, then e w 0.004; if hadronic processes 
contribute to the gamma-ray luminosity, this is an upper 
limit. 

The synchrotron luminosity also depends on the mag- 
netic field strength, and we consider two physical scenar- 
ios for its strength and evolution. In the first case, we as- 
sume that the magnetic field in the forward-shock region 
is simply the magnetic fie ld from the Mira wind, com- 
pressed. As suggested bv lAbdo et al.l ()2010D . the mag- 
netic field energy density in the wind is expected to be 
in equipartition with the thermal energy density of the 
wind. We assume that the relevant conditions pertain to 
the wind before it was heated and ionized by the nova 
outburst, so we set the wind temperature to Tu,o — 700 
K as in Abdo et al. The magnetic field strength is then 
simply a function of distance from the Mira giant, and is 
described as: 

-^comp — v32 Tt-ncsukT^Q. (13) 

Another possibility is that the magnetic field is am- 
plified in the post-shock region to be in equipartition 
with the relativistic electrons. This simplification is com- 
monl y used to interpret th e radio emission from SNe 
Ce.g.. iChevalier fc Franssonl [2006). and was recently uti- 
hzed bv lChugail ^2mW for outbursts in symbiotic stars. 
In this case: 

Samp = \/8TTeiJ,mH TlCSM v"^ . (14) 

Here, B depends on both the distance from the Mira 
giant and the blastwave velocity. 

The synchrotron flux density of a spherical body of 
radius R at distance D (where both R and D are in 
units of cm) can be found by Equation 1 of iChevalieii 
(IT998I) to be: 



= 5.64 X 10 



-22 



1/2 „5/2, 



1 — exp 



C6 

-(p+4)/2- 



{R/DY mJy 



(15) 



Here, C5 and cg are constants that depend on p as tabu- 
lated in [Pacholczvk (,197Qi ). and as listed in Table |3] for 
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Figure 10. Models for synchrotron emission from V407 Cyg, overplotted as dashed lines on our measured multi-frequency light curve 
(solid lines and diamonds). The two panels represent di ffere nt assumptions about the magnetic field strength i n the forward-shocked region: 
the left panel assumes -Bcomp as described in Equation 1131 while the right panel assume Bamp from Equation 1141 



convenience. The observing frequency is v in units of Hz, 
and 1^1 is defined as: 



i^i = 1.25 X 10 
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2/(p+4) 



^(p+2)/(p+4) 

(16) 

We assume that the synchrotron emission is contained 
in a spherical shell between the forward shock and the 
contact discontinuity, typicall y with a width ^ 20% of the 
contact discontinuity's radius ()Chevalieiill982[ ). In Equa- 
tion [161 / is the volume filling factor of this emission, so 
/ = 0.42 for a width of 20%. 

All the above equations assume spherical symmetry, 
but we know that V407 Cyg is an inherently asymmetric 
case. Because synchrotron luminosity is sensitive to both 
shock velocity and circumstellar density, it is important 
that we take into account this asymmetry in predicting 
the expected synchrotron si gnal from V407 Cy g. We 
again use the toy model of iNelson et al.l (j2012t ) to es- 
timate the evolution of the blastwave in this asymmet- 
ric medium (see Section I5.1.2p . We calculate the syn- 
chrotron flux density emanating from each polar-angle 
ring by using Equation [15] and scaling the flux density 
by the fraction of a spherical shell's volume that is en- 
closed in that ring (/vol! Equation [3]) . The fiux densities 
from all rings are summed together to give the total syn- 
chrotron flux at a given time, producing the model light 
curve plotted in Figure [TUl 

The light curves in Figure [TOj show the evolution of 
this synchrotron model, assuming Mej = 10^^ Mq, 



= 2 X 10-^ Mq yr-i, = 10 km s'^, e = 0.004, 
a = 17 AU, an initial shock velocity of Vg = 3, 200 km 
s^^, and a distance of 3.0 kpc. The left panel assumes a 
compressed magnetic field (Equation [TS]). while the right 
panel assumes an amplified magnetic field in equiparti- 
tion with the relativistic electrons (Equation (Til) . The 
synchrotron emission in both models is optically thin, 
with spectral index a = —0.65, so the lower frequen- 
cies are brighter. The synchrotron flux density decreases 
with time because the blastwave is decelerating and the 
density of the circumstellar material is declining in all 
directions that are significantly contributing to the syn- 
chrotron emission. The predicted synchrotron emission 
in the -Bcomp case is much less luminous than in the Bamp 
case, because the magnetic field in the forward shock is 
weaker. 

Both models predict bright synchrotron emission that 
would have been easily detectable with the VLA in the 
first months after outburst, with flux de nsities >> 1 mJy 
(a re sult confirmed by the models of iMartin fc DubusI 
l2012t ). However, in reality at early times, when the 
synchrotron luminosity was highest, the synchrotron- 
emitting forward shock suffered free- free absorption from 
the ionized Mira wind (compare the black solid line in 
Figure [7| with the radii of t = 1 surfaces). The large 
radius of the absorbing screen means synchrotron emis- 
sion is only visible at a given frequency once the wind 
becomes optically thin at that frequency. A compari- 
son between synchrotron models and observations is most 
simply carried out on Day ^^470, when the observed emis- 
sion from V407 Cyg is consistent with optically thin ther- 
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mal radiation at all observed frequencies. In the -Bcomp 
model, synchrotron emission constitutes a trivial fraction 
of the total flux density on Day 470, even at 1.4 GHz. In 
the -Bamp model, on the other hand, the predicted syn- 
chrotron luminosities at lower frequencies (< 8 GHz) are 
significantly higher than those observed on Day 470, and 
we might have expected to see a clear signature of syn- 
chrotron emission in the spectral evolution. We stress 
that no such synchrotron contribution is observed in our 
radio spectra; spectral indices measured between Days 
359-669 are a w 0, as expected for optically-thin thermal 
emission {a = —0.1) and inconsistent with expectations 
for synchrotron emission (a « —0.65). We may be able 
to constrain the magnetic field in the forward shock with 
our radio light curve. 

Of course, our model of the synchrotron emission de- 
pends quite sensitively on a number of assumed physical 
parameters. Although V407 Cyg is a better-constrained 
system than many, thanks to the community's wealth of 
multi- wavelength data, there is still room to vary predic- 
tions of the synchrotron luminosity by orders of magni- 
tude. For example, a decrease of the ejecta mass from 
Mej = 10"^ Mq to Mcj = 10""^ Mq means that the 
blast wave decelerates much faster, and decreases the 1.4 
GHz fiux density predicted by the i?amp model on Day 
470 from 11.5 mJy to 0.2 mJy. Because synchrotron lu- 
minosity depends so sensitively on blast wave velocity, 
we are likely too simplistic in assuming that the ejecta 
are expelled in a thin shell with a single velocity; the 
synchrotron model might be improved with more realis- 
tic radial distributions of ejecta mass and velocity. The 
Samp model depends on shock velocity more sensitively 
than the -Bcomp model; decreasing M^j to 10~ Mq in the 
Scomp case only decreases the 1.4 GHz fiux density from 
0.55 mJy to 0.11 mJy. In the Samp case, the energy den- 
sities in both the magnetic field and relativistic electrons 
depend on the shock velocity, while in the i?comp case, 
the magnetic field is independent of shock velocity. 

The synchrotron luminosity also depends sensitively 
on assumptions about the efficiencies of particle accel- 
eration and B field amplification in the shock. Al- 
though V407 Cyg is unusual because e has been directly 
constrained from the gamma-ray observations, this con- 
straint remains a rather uncertain upper limit. Ad- 
ditional uncertainties are involved in our assumptions 
about the evolution of e with time and the degree of 
equipartition between relativistic electrons and magnetic 
field. In the i?amp case, if e = 0.002 instead of e = 0.004, 
the 1.4 GHz fiux density on Day 470 would have been 
3.2 mJy rather than 11.5 mJy. Assuming the magnetic 
field and relativistic electrons are in equipartition, the 
magnetic field in the shock on Day 470 must have been 
< 0.009 G to produce less than the observed 1.4 GHz 
flux density. 

While the predictions for synchrotron luminosity re- 
main uncertain, reasonable parameters imply that signif- 
icant synchrotron emission was produced at early times, 
but it was hidden behind the r = 1 surface of the ionized 
Mira wind for most of the evolution of the nova. Despite 
rich evidence for particle acceleration in V407 Cyg, it is 
not surprising that we did not detect synchrotron emis- 
sion in our VLA data. In contrast, if the thermal radio 
emission had been dominated by the nova ejecta, the syn- 



chrotron emission from the forward shock would not have 
been absorbed, and would have been detectable — yet an- 
other reason for dismissing the thermal ejecta scenario 
investigated in Section |4] as the source of radio emission 
in V407 Cyg. 

7. HISTORICAL RADIO OBSERVATIONS OF V407 CYG 

Although only radio non-detections of V407 Cyg have 
been been reported in the literature for the time period 
before the 2010 nova outburst (Table |4]), we used unpub- 
lished archival data from the VLA to determine that the 
source did in fact produce detectable radio emission in 
1993, near the begi nning of a period of increased activity 
of the WD (Koloti lov et al.| [l998'). 

The archival VLA data, obtained on 1993 May 14 as 
part of program AI48 (PI R. Ivison), were acquired in 
BnC conflguration at X band with 100 MHz bandwidth; 
they were calibrated using 3C286 and J2007-I-4029. With 
15 minutes on source, the 1993 data yield a fiux density 
for V407 Cyg of 1.18 ± 0.07 mJy at 8.4 GHz. This de- 
tection is in contr ast to the non-detectio n obtained four 
years earlier, when'Sea auist et al.l (|1993[ ) found an upper 
limit on the fiux density of 0.06 mJy at 8.4 GHz. Thus, 
the strength of the radio emission from V407 Cyg varied 
by more than a factor of 20 even before the 2010 nova 
eruption. 

The radio luminosity of V407 Cyg probably increased 
around 1993 because of a fiare of ionizing radiation from 
the WD. From spectroscopic monitoring of V407 Cyg 
over the last two decades, it is known that the luminosi- 
ties of optical emission lines vary by factors of ~5 over 
just a few months. On 1993 May 29, just two weeks after 
our radio detection, an optical spectrum of V407 Cyg did 
not show high-excitation signatures of strong symbiotic 
activity (although it did feat ure Balmer emission, w hich 
had strengthened since 1991: iKolotilov et al.llTgM ). By 
August 1994, emission lines of high ionization states (e.g., 
[O III]AA4959,5007, He IIA4686) were visible in the op- 
tical spectrum for the first time si nce modern observ- 
ing began of this so urce in 1984 (jMunari et al.l [1991 
IKolotilov et al.iri998f ). We hypothesize that this active 
period of V407 Cyg may have begun in early 1993 and 
gradually increased in power. The radio appears to pre- 
cede the strongly symbiotic optical spectrum, but note 
that in spring 1993, the Mira giant was near the maxi- 
mum of its pulsation cycle, and the giant's light would 
have diluted the optical emission line spectrum. In Au- 
gust 1994, the giant is approaching a pulsational mini- 
mum, and the relatively faint higher-excitation emission 
lines would have been more easily visible. 

Several years late r, in 1998, an ou tburst was observed 
in V407 Cyg (Shug arov et al.ll2007[ ). Unfortunately, no 
radio light curve exists for the interesting period be- 
tween 1993 and 1998; V407 Cyg was not observed again 
with the V LA (except with a shallow NVSS exposure 
at 1.4 GHz: iCondon et al.|[T99l Table HI) until the 2010 
nova. In addition, no information on the spectral index 
exists for the 1993 detection of V407 Cyg, so the ra- 
dio emission mechanism can not be confirmed. In other 
symbiotic binaries, dramatic increases in radio fiux have 
been associated with collimated mass ejection, which 
can produce both thermal and sync hrotron radio emis- 
sion (e.g., CH Cyg, R Aqr, and M ira; 'Taylor et al."1986l 
iCrocker et all 120011: iNichols et al. ,2007; .Meaburn et ail 
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I2009D . While there is, to date, no evidence for last 
coUimated outflows in optical spectroscopic monitoring, 
V407 Cyg remains a promising candidate site for the pro- 
duction ol synchrotron emission. The dramatic radio 
variability of V407 Cyg calls for long-term monitoring 
of the radio flux density and spectral index, in coopera- 
tion with synoptic optical spectroscopy, in order to better 
understand symbiotic activity in this source. 

8. DISCUSSION AND CONCLUSIONS 

Our VLA light curve, spanning frequencies 1.4-45 GHz 
and covering 17-770 days after outburst, shows a radio 
evolution that is unusually luminous and slowly-evolving, 
compared with other novae. Given the postulated low 
ejecta mass and lack of detected synchrotron emission, 
the most plausible source of radio emission from the nova 
event in V407 Cyg is the ionized wind from the Mira gi- 
ant. The radio emission from other symbiotic binaries is 
explained with a similar mechanism, attributing free- free 
emission to a wind t hat has been ionized by a hot WD 
(jSeaquist et al.lflQS^ . However, the radio emission from 
V407 Cyg is unique due to the relatively rapid rise and 
decline of this radio emission. To achieve a sudden rise 
in the radio light curve requires an abrupt and strong 
increase in the ionizing flux, and the dramatic fading of 
the radio light curve requires a mechanism for disrupt- 
ing or heating the wind on relatively short timescales. 
We hypothesize that the radio light curve declines as 
the wind slowly recombines on the outside and is shock- 
heated from within by the nova blast. The radio light 
curve provides an estimate of the mass-loss rate of the 
Mira, M„ w 2 x 10"^ Mq yr'^ (for Vy, = 10 km s"^ and 
d — 3 kpc) . We place a relatively hard lower limit on the 
mass-loss rate, assuming the wind is completely ionized, 
of M„ > 7 X 10"'^ Mq yi-\ 

8.1. Different Observational Tracers Imply Discrepant 
Mira Mass Loss Rates in V407 Cyg 

A few difficulties remain for our claim that « 
10~^ M0 yr~^ in V407 Cyg. In a the Suzaku X-ray spec- 
trum of V407 Cyg observed on Day 30, non-equilibrium 
ionization lines are observed, implying an ionizat ion age 
in th e forward shock of 2 x 10^° cm~^ s (Nelso n~et al.l 
l2012t ). This relatively small ionization age translates to 
an upper limit on the density of circumstellar material 
at T-RG ~ 10^^ cm of ncsM ^ 10^ cm~^ — or a mass-loss 
rate ^ 4 x 10~^ yr~^ (an order of magnitude 

lower than determined from the radio; both assume 
Utu = 10 km s~^). In addition, studies of the hard X-ray 
light curve can successfully model the data as shocked 
Mira wind, but M,,, » 10~^ - 10"''' M(7) yr~ ^ is required 
(jOrlando fc Drakel[20Tl iNelson et al]l2012l ). If the wind 
is too dense, the emission measure of the forward-shocked 
material will be too high and the ejecta will decelerate 
too rapidly, making the hard X- ray light curve pea k ear- 
lier and brighter than observed ()Mukai et al.ll2012l ). 

A similar effect can be seen in the maximum expansion 
velocity as a function of time, measured from the wings 
of Ha emission line profiles (Sh ore et al.ll20lj] ). If the ra- 
tio of the wind density to the ejecta mass is too high, the 
ejecta will decelerate faster than observed. We compare 
the redshifted measurements from Shore et al.'s Figure 
19 with our simple model of the asymmetric blastwave 



(|Nelson et al.llMl) . The nova material traveling away 
from the giant companion maintains the highest velocity, 
so we compare this direction with the Shore et al. data. 
We find that the observed velocity evolution can be re- 
produced if, for an ejecta mass of Mcj — 10~^ Mq, the 
wind mass loss rate is M^, « 3 x 10~* yr~^ (again, 
an order of magnitude lower than our radio-determined 
mass- loss rate). 

On the other hand, the infrared spectral energy dis- 
tribution of V407 Cyg implies « 8 x lO""^ M© yr"! 
(jKolotilov et al.iri998f ). in relatively good agreement with 
our radio estimate. 

All estimates of in V407 Cyg assume that the 
Mira wind is spherically symmetric, which is likely 
an over-simplification. Simulations show that bi- 
nary interactions lead to a concentration of circumbi- 
nary material in the orbital plan e ([Waldcr ct al. 200^ 
IMohamed k, Podsiadlowskil 1201 1[ ). lOrlando fc Drakg 
(I2012D find that such a density enhancement is critical 
for matching the X-ray light curve of V407 Cyg with a 
small recurrent-nova-like ejecta mass. Radio imaging of 
V407 Cyg with MERLIN and VLA over the weeks to 
years following the 2010 nova sho w structure in the cir- 
cumbinary medium (|O'Brienll2010l O'Brien et al. 2012, in 
preparation, Mioduszewski et al. 2012, in preparation). 
However, the fact that the observed radio luminosity is 
primarily dependent on the area of r = 1 surfaces implies 
that we can not significantly decrease the radio estimate 
of Myj by summoning large-scale structure in the wind. 
For example, consider the case of a wind funneling most 
of its mass into a disk-like structure; if the disk is viewed 
edge-on, the r = 1 surface will be relatively small and the 
radio emission will be faint (compared with a spherical 
wind). The material in the disk will be dense, implying 
that, if we view it face on, it will only be ionized out to 
relatively small radius (again producing small t ~ 1 sur- 
faces) and will recombine quickly. For these reasons, the 
most effective way to produce large radio luminosities is 
to distribute the wind roughly spherically. 

Small-scale "clumps" in the Mira wind can boost the 
emission measure of the wind and lead to over-estimates 
of from radio data CA bbott et allll981l: iNugis et al.1 
|1998() . Small-scale clumping, on the other hand, should 
not affect the IR, X-ray, or Ha dynamical estimates of 
Mw SiO mas ers are detected in the env elope surround- 
ing V407 Cyg peguchi et al.ll2005l[20TTl ). a p henomenon 
that likely requires an inhomogeneous wind (jGrav et al.l 
119981) . In addition, the X-ray and radio light curves of the 
2010 nova show rapid variations (Figure [T| INelson et al.l 
|2012| ) , which might be produced by structure in the wind. 
This variety of multi-wavelength diagnostics probe dif- 
ferent radii in the wind (^ 10^^ cm from the Mira for 
masers; ~ 10^^ cm for the X-ray light curve; ~ 10^^ — 10^^ 
cm for the radio light curve), so we see evidence for in- 
homogeneities at many radii in the wind of V407 Cyg. 
Therefore, it is possible that our radio measurement of 
Mw is an over-estimate, and it might be brought into 
better alignment with the X-ray and Ha dynamical es- 
timates by accounting for clumping. However, if the fill- 
ing factor of the wind is significantly less than unity, the 
ionization state of the circumbinary material could dif- 
fer significantly from what we predict in Section [Sj as it 
would recombine much more quickly. 
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Observed mass loss rates of Mira giants span M^^ = 
few X 10~* — few x 10""* Mq yr~^, with an average 
(M„,) = 6 X 10-'^ Mq yr-i (IKnapp fc Morris! [19851 al- 
though this sample is not complete for the lowest M^, 
stars). The Mira giant in V407 Cyg has a very long pe- 
riod and is oxyge n- rich, so it is more typical of OH/IR 
stars than Miras (jMunari "eTall [199(11) . OH/IR sources 
typic ally show very de nse winds with ~ 10"*^ M© yr~^ 
(e.g., iBaud k Habh^[T98a) . but V407 Cyg is an odd 
case because dust formation in the wind is suppressed 
by the hot WD, and therefore the wind-d riving mecha- 
nism may differ from typical OH/IR stars (jMunari et al.l 
[l99l . Still, we might expect V407 Cyg to occupy the 
upper end of the Knapp & Morris range, rather that the 
lower end — more consistent with the radio and IR mass- 
loss rates than those implied by the X-ray data and the 
deceleration of the nova ejecta. 

We therefore postulate that the mass-loss rate in 
V407 Cyg is ~ 10~^ Mq yr~^. One possible solution for 
reconciling the radio-derived M^, with the lower mass- 
loss rates implied by X-ray non-equilibrium ionization. 
X-ray emission measure, and shock deceleration is a cir- 
cumbinary medium that is concentrated in the orbital 
plane at small radii (rwD ^ 10"'^^ cm) and becomes more 
spherically distributed at larger radii (rwD ^ 10^^ cm). 
Such a non-spherical wind is in rough kee ping with the 
simula tions of circumbinary material by iWalder et al.l 
(|2008D . which show an over-density in the equatorial 
plane of a factor of ~ 10 — 100 for r-yvD ^ a, and a 
flare in the outer circumbinary disk so that the equa- 
torial over-density at larger radius is only a factor of 
^2 — 3. However, this scenario would require significant 
fine tuning to come into line with the multi-wavelength 
data on V407 Cyg; the X-ray ionization timescale implies 
a low density at ^ 10^^ cm while the 25 GHz light curve 
measures ~ 10~^ Mq yr~^ at ~ 10^^ cm. A sudden 
radial transition would need to occur in the distribution 
of circumbinary material, from a "disky" concentration 
in the equatorial plane to a near-spherical distribution 
at a radius of ~ 10^^ cm to reconcile these data. The 
X-ray and radio data may be brought into better align- 
ment if the binary separation is significantly larger than 
17 AU (in wh ich case, the estima ted orbital period of 
43 years from [Munari et all 119901 would be incorrect). 
A larger orbital separation would decrease the density 
of circumbinary material away from the red giant while 
still allowing for a relatively high M^,. Our model for the 
radio light curve would be largely unchanged, but it re- 
mains unclear if the X-ray light curve could be explained 
with such a wide binary separation ('Orlan do fc Draki 
120121 kept a fixed at 15.5 AU in their detailed modeling 
of the X-ray light curve). 

8.2. Implications for Radio Emission from Novae 

The physical process that we claim is responsi- 
ble for the radio light curve of V407 Cyg — namely, 
the time-variable ionization state of the circumbinary 
environment — has never before been recognized as the 
primary source of radio emission in a nova outburst. 
It provides an interesting test case of how dense cir- 
cumstellar material can affect the evolution of a radio 
transient. Similar outbursts of thermal radiation, pro- 
duced by the flash ionization of circumbinary material. 



might be visible in other novae if significant masses of 
circumbinary mate r ial ar e prevalent, as suggested by 
IWilhams fc ([20Tol) . Of course, this emission — 

particularly if it is optically thin — might be superim- 
posed on, or drowned out by, other emission compo- 
nents, like the thermal radiation from the expanding 
nova ejecta or synchrotron emission from the nova blast. 
The most similar known nova to V407 Cyg, RS Oph, 
was surrounded by a lower density wind; it did not form 
an absorbing screen following the nova flash, and there- 
fore, synchrotr on emission was direc t ly detectable from 
the nova blast CHicllming ct al."1986': jTavlor et al.lll989l : 
lO'Brien et all r200& .Rupen et al, 2QM)- The recurrent 
nova T Pyx is surroun ded by an Ha-|-[N II] nebula of 
mass ^ few x 10"^ Mq (|Shara et al.lll997t iSchaefer et al.l 
l2010f ) which might have been flash-ionized by the 2010 
nova outburst. However, scaling from the Hq;+[N II] lu- 
minosity of this nebula in the months after outburst, we 
find that the nebula produced a negligible <1 /iJy of ther- 
mal radio emission (Nelson et al. 2012, in preparation). 
Therefore, V407 Cyg remains the only nova to date that 
shows a radio light curve dominated by the temperature 
and ionization structure of its circumbinary medium. 

8.3. Implications for the Progenitors of SNe la 

The small nova ejecta mass postulated for V407 Cyg 
implies that the WD is massive (>1.25 Mq; 
iMik olaiewskal [2010l ). making this system an intriguing 
progenitor candidate for a SN la. A symbiotic progenitor 
like V407 Cyg was conclusi vely ruled out in the recent 
case of SN2011fe in MlOl (ILi et al.ll20lH lHoresh_eLaD 



I2012t iChomiuk et"al1 l2012t iMargutti et a l. 2012), and 
has also been ruled out for an additional '^15 SNe with 
deep VLA limits (Chomiuk et al. 2012, in preparation). 
However, the sustained elusiveness of SN la progenitors 
and subtle diversity of SN la explosions has led the com- 
munity to consider the possibili ty that SN la pro genitor 
systems may be diverse (e.g., iFolev et all l2012| ). For 
example, while most SNe la show no evidence f or hydro- 
gen i n their spectra down to very deep limits (jLeonardl 
[200l . a very small handful of SNe la develop luminous 
Ha emission lines in their spectra, indica tive of strong 
inter a ction with circum stellar material (|Hamuv et al.l 
120031 : iDildav et al.l [20121 ). In addition, while some SNe 
la show time-variable blue-shifted Na I D absorption 
lines, argued to be evidence of circumstellar material, 
other s do not (jPatat et all 120071 . 12011bl: [Sternberg erahl 

MM- 

Perhaps a fraction of SNe la have symbiotic progen- 
itors, with varying delay times between mass loss and 
SN explosion. For those few WDs that explode into a 
gas-rich medium like present-day V407 Cyg, the explo- 
sion may resemble the Ha-bright SNe la. More com- 
monly, SN la-progenitor WDs are surrounded by low- 
density environments; their immediate environment may 
have been s wept clean by nova blasts that swept up the 
giant wind ([Patat et al.ll2011at iWood-Vasev fc Sokoloskil 
l200(]f ). However, such a scenario does not describe 
present day V407 Cyg because, using the classical re- 
lation between radio fiux density and Mu, (Equation 
[TU|) . we derive very similar mass- loss rates for V407 Cyg 
for frequencies spanning the range 1.4-25 GHz, which 
are sensitive to the wind density at very different radii 
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lO^'^ — 10"'^^ cm). These observations imply that there 
is no dramatic radial structure in the wind for radii 
> 10^^ cm, and the Mira envelope follows a smooth 
ncSM fx profile. 

Finally, the circumbinary medium around a symbiotic 
system might have time to expand away and merge with 
the ambient interstellar medium, if there is a sign i ficant 
delay between accreti on and explosion pusthami 120111 : 
iDi Stefano et al.ll201lD . V407 Cyg could constitute such 
a progenitor in ^ 10^ years time. In this case, any blue- 
shifted Na I D absorption lines present in the resulting 
SN la spectrum would be the product of the Mira enve- 
lope slowly drifting away from its parent binary. 
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VLA Observations of V407 Cyg 
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1.8 


8.23 ±0.51 


19.0 


14.87 ± 1.49 


25.0 


13.63 ± 1.37 


1.4 


7.92 ± 0.52 


1.8 


8.19 ±0.46 


19.0 


13.26 ± 1.33 


25.0 


13.30 ±1.33 


1.4 


5.49 ± 0.38 


1.8 


6.39 ±0.40 


19.0 


11.87 ± 1.19 


25.0 


11.20 ±1.12 


1.4 


5.61 ± 0.37 


1.8 


6.34 ±0.36 


19.0 


10.01 ± 1.00 


25.0 


9.51 ±0.95 


1.4 


5.58 ± 0.31 


1.8 


5.53 ±0.31 


19.0 


9.86 ± 0.99 


25.0 


10.08 ± 1.01 


1.3 


4.22 ± 0.27 


1.8 


4.64 ±0.25 


20.1 


4.65 ± 0.47 


25.6 


4.74 ±0.48 


1.3 


3.60 ±0.20 


1.8 


4.07 ±0.21 


20.1 


2.54 ±0.26 


25.6 


2.62 ±0.27 


1.3 


3.64 ±0.56 


1.8 


3.86 ±0.35 


20.1 


3.67 ±0.37 


25.6 


3.83 ±0.39 


1.3 


4.69 ± 1.00 


1.8 


3.32 ±0.55 


20.1 


3.14 ±0.32 


25.6 


2.99 ±0.30 


1.3 


2.94 ±0.33 


1.8 


3.35 ±0.29 


20.1 


2.91 ± 0.29 


25.6 


2.97 ±0.30 



5.0 


« QQ -L n 
u.oy HZ Kj.oo 






5.0 


7.03 ±0.37 






4.8 


f; CO _L ri 07 


7.0 


0.0 ± in Kj.^o 






49.0 


44.98 ± 22.55 


4.8 


Q QQ _L f) 


7.0 


11 fts -u (1 fin 


4.8 


1 f) «o _|_ f) c;c; 

J-W.UO _1_ \J.iJiJ 


7.0 


J- z . 1 _1_ VJ.UIJ 


41.0 


59.89 ± 12.01 


45.0 


65.63 ± 13.18 


4.8 


12.36 ±0.65 


7.0 


14.67 ±0.74 


4.8 


12.79 ±0.66 


7.0 


15.59 ±0.79 


4.8 


13.69 ±0.71 


7.0 


15.72 ±0.79 


4.8 


13.96 ±0.70 


7.0 


16.07 ±0.81 


4.8 


13.67 ±0.69 


7.0 


16.88 ±0.85 


4.8 


13.14 ±0.67 


7.0 


16.10 ±0.81 


4.8 


14.16 ±0.71 


7.0 


17.18 ±0.86 


4.8 


15.02 ±0.76 


7.9 


19.18 ±0.96 


41.0 


51.34 ± 10.30 


45.0 


50.66 ± 10.16 


4.8 


15.48 ± 0.79 


7.9 


19.88 ± 1.00 


41.0 


44.90 ± 9.00 


45.0 


46.31 ±9.28 


4.8 


17.40 ±0.88 


7.9 


21.71 ± 1.09 


41.0 


39.97 ±8.01 


45.0 


38.40 ± 7.71 


4.8 


17.47 ±0.89 


7.9 


20.62 ± 1.04 


41.0 


37.05 ± 7.43 


45.0 


35.88 ±7.20 


4.8 


17.26 ±0.87 


7.9 


20.19 ± 1.01 


4.8 


17.78 ±0.90 


7.9 


19.02 ±0.95 


4.8 


17.45 ±0.88 


7.9 


18.20 ±0.93 


41.0 


25.94 ±5.23 


45.0 


28.13 ±5.68 


4.8 


12.93 ±0.65 


7.9 


14.57 ±0.73 


41.0 


18.67 ±3.83 


45.0 


18.61 ±3.85 


4.8 


11.81 ±0.60 


7.9 


12.92 ±0.65 


41.0 


13.39 ±2.52 


45.0 


13.51 ±2.75 


4.8 


10.79 ±0.55 


7.9 


11.31 ±0.57 


41.0 


14.58 ±2.97 


45.0 


15.12 ±3.08 


4.8 


9.98 ±0.51 


7.9 


10.16 ±0.51 


41.0 


9.66 ± 2.11 


45.0 


9.09 ± 1.98 


4.8 


9.17 ±0.46 


7.9 


9.10 ±0.46 


41.0 


14.19 ±2.91 


45.0 


12.75 ±2.64 


5.0 


5.00 ±0.25 


7.4 


4.81 ±0.24 


41.0 


5.30 ± 1.08 


45.0 


5.45 ± 1.11 


5.0 


3.27 ±0.17 


7.4 


2.75 ±0.14 


41.0 


4.25 ±0.87 


45.0 


4.45 ±0.92 


5.0 


3.54 ±0.18 


7.4 


3.49 ±0.18 


41.0 


4.32 ±0.87 


45.0 


4.47 ±0.91 


5.0 


3.21 ±0.17 


7.4 


3.15 ±0.16 


41.0 


3.96 ±0.83 


45.0 


5.23 ± 1.11 


5.0 


2.91 ±0.15 


7.4 


2.80 ±0.14 


41.0 


2.88 ±0.60 


45.0 


2.66 ±0.59 



Time since optical discovery of the nova event (2010 March 10.8: IMunari et al.ll201ll ). 



Table 2 

Spectral Indices for V407 Cyg 







1-45 GHz 


1-9 GH 


z 


19-45 GHz 


Epoch 




a 


Cdof ~) 


o 


X^(dof) 


ol 


^P■ f dof 1 


1 b 

-i 


11.7 


1.03 ± 0.18 
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0.98 ±0.26 


3.45(2) 


o 
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i.y.u 


1.34 ±0.29 





l.o4 ± U.zy 


U 






Q 

o 


/o.y 


0.75 ± 0.02 


0.28(5) 


n ^^c J— o ni 
U.Do ± U.Ui 
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U.Ui^z ) 


0.73 ±0.10 


0.06(1) 


4 
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0.69 ± 0.01 


0.06(2) 


n '7 A J— O 1 o 

0. /4 ± u.iy 





0.58 ±0.52 










0.76 ±0.02 


0.50(6) 


n 7/1 J— o nc; 
U. / 4 ± U.UO 


U.oyi^z ) 


0.94 ±0.08 


0.13(2) 


D 


oo.y 


0.77 ± 0.05 


2.06(4) 


o QO ^ o no 




0.68 ±0.52 





7 


38.9 


0.70 ± 0.03 


1.17(4) 


o '7R _L n oc; 


0.91(2) 


0.60 ±0.52 





8 


42.9 


0.61 ± 0.05 


1.42(4) 


o _L n oc; 
U.6d ± O.Uo 


1.29(2) 


0.17 ±0.52 





n 

y 


40. y 


0.62 ± 0.03 


0.79(4) 


o P.^ _L n oc; 
U.DO it U.UO 


1. i4(^Z j 


0.39 ±0.51 







/in Q 


0.64 ±0.01 


0.17(4) 


o c\A _L n oo 
U.04 it U.UZ 


n o'7/'o\ 
U.z / l^zj 


0.46 ±0.52 





1 1 




0.61 ± 0.02 


0.40(4) 


o fil 4- n 0/1 
U.Ui It U.U4 


U. 


0.65 ±0.52 





1 o 


c;'7 Q 
o / .6 


0.55 ± 0.03 


1.05(4) 


O _L O oo 

U.Oi it U.Uz 


n 1 A/'o\ 
U. iul^z ) 


0.22 ±0.52 







Oo 


0.63 ± 0.04 


3.01(6) 


o '7/1 _L o oc; 
U. / 4 it U.UO 


1. r o(^z ) 


0.72 ±0.16 


0.54(2) 


1 A 


DO. 6 


0.58 ± 0.04 


3.01(6) 


O 70 _L O Ot; 

U. / U it U.UO 


1.4 / (^Zj 


0.56 ±0.17 


0.58(2) 




r o.o 


0.56 ±0.06 


6.56(6) 


O 7/1 _L O 07 

U. / 4 it U.U / 


o.z4(^z j 


0.39 ±0.12 


0.30(2) 


16 


79.8 


0.55 ± 0.06 


7.05(6) 


O 7Q _L O oo 

U. / o ± U.U8 


4.33(2) 


0.37 ±0.10 


0.20(2) 


17 


96.7 


0.74 ± 0.09 


5.38(2) 


o 7/1 _L n on 
U. / 4 ± U.U9 


5.38(2) 






18 


101.8 


0.68 ±0.12 


8.30(2) 


0.68 ± 0.12 


8.30(2) 






19 


102.7 


0.47 ± 0.07 


9.26(6) 


0.67 ±0.11 


8.31(2) 


0.20 ±0.10 


0.21(2) 


20 


138.6 


0.27 ± 0.05 


4.73(6) 


0.43 ± 0.04 


1.17(2) 


0.27 ±0.21 


0.80(2) 


21 


153.6 


0.19 ±0.04 


2.38(6) 


0.30 ±0.03 


0.54(2) 


-0.03 ±0.05 


0.05(2) 


22 


168.5 


0.28 ±0.05 


4.88(6) 


0.43 ±0.07 


2.89(2) 


0.27±0.17 


0.54(2) 


23 


183.5 


0.20 ±0.05 


5.41(6) 


0.35 ±0.06 


2.62(2) 


-0.09 ±0.05 


0.04(2) 


24 


203.4 


0.25 ± 0.04 


2.82(6) 


0.33 ±0.07 


3.85(2) 


0.36 ±0.12 


0.31(2) 


25 


359.6 


0.04 ± 0.02 


0.56(6) 


0.07 ±0.03 


0.82(2) 


0.19 ±0.03 


0.02(2) 


26 


469.6 


-0.13 ±0.05 


3.13(6) 


-0.17 ±0.07 


3.86(2) 


0.71 ±0.19 


0.59(2) 


27 


603.3 


0.02 ± 0.03 


0.45(6) 


-0.05 ±0.02 


0.10(2) 


0.24 ±0.02 


0.01(2) 


28 


669.0 


-0.01 ±0.05 


1.15(6) 


-0.10 ±0.07 


0.62(2) 


0.47 ±0.26 


1.04(2) 


29 


769.6 


-0.03 ±0.02 


0.40(6) 


-0.07 ±0.05 


0.61(2) 


-0.07 ±0.24 


0.06(2) 



Time since optical V-band maximum (2010 March 10.8: IMunari et aI.ll20Tll '). 
Prom 2010 March 22/23, as documented in lNestoras et ai.l ll2010f) . 



Table 3 

Constants for Calculating Synchrotron 
Emission (Pacholczyk 1970) 



V C5 C6 



2.0 7.52 X 10"^"' 7.97 x 10" 
2.5 9.68 X 10-24 8.10 x 10" 
3.0 1.37 X 10-23 8.61 X IQ- 



Table 4 

Historical Radio Observations of V407 Cyg 



Date 


Frequency 




Facility 


Reference 


(UT) 


(GHz) 


(mJy) 






1982 Jul 13 


4.9 


< 0.96 


VLA 


Seaauist et al. (1984) 


1988 Oct 16 


5.0 


< 1.72 


BBP 


Ivison et aL (1991) 


1989 Dec 22 


8.4 


< 0.06 


VLA 


Seaquist ct al. (1993J 


1992 Oct 10 


270 


< 114 


JCMT 


Ivison et al. (1995) 


1993 May 14 


8.4 


1.18 ±0.07 


VLA 


This Work 


1996 Aug 18 


1.4 


< 2.77 


VLA 


Condon et al. C 1998) 



" The Jodrell Bank Lowell-Mark 2 Broad Band Interferometer. 



